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ABSTRACT 



Perfect bubbles like that surrounding the galactic Hii region RCW 120 (jPeharveng et al 

20091) have been interpreted as proof of concept for the collect and collapse (C&C) mechanism 
of triggered star formation. The cold, dusty clumps surrounding RCW 120 are aligned along 
an almost spherical shell. It has been inferred that these massive clumps, which sometimes 
harbour young stellar objects, have been formed via the fragmentation of the dense, swept-up 
shell. In order to better understand the triggering mechanisms at work in shells like RCW 120, 
we perform high-resolution, three dimensional SPH simulations of HII regions expanding into 
fractal molecular clouds. In a second step, we use RADMC-3D to compute the synthetic dust 
continuum emission from our simulations, in order to compare them with observations of RCW 
120 made with APEX-LABOCA at 870 ^m. We show that a distribution of clumps similar to 
the one seen in RCW 120 can readily be explained by a non-uniform underlying molecular cloud 
structure. Hence, a shell-like clump configuration around an HII region does not necessarily 
support the C&C scenario, but rather reflects the pre-existing, non-uniform density distribution 
of the molecular cloud into which the HII region expands. 
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1. Introduction 

The possibility of triggering star formation by 
ionizing feedback from young, massive stars has 
been explored for several decades. From a theo- 
retical point of view, two main triggering mecha- 
nisms have been suggested: Collect and Collapse 
(C&C), and Radiation Driven Implosion (RDI). 

The CfcC mechanism was first analyzed by 
Elmegreen & Ladal (jl977() . In this mechanism the 
expanding HII region sweeps up a layer of cold 
gas and dust beyond the ionization front (e.g. 
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Dale et al] l2007l) . and this shell eventually be- 
comes gravitationally unstable due to the growth 
of pe rturbations along its surface jElniegreer 



Wunsch et al 



1994 IWhitworth et all Il994al: iDale et alj 
l2010[) 



One argument in favour 
of the C&C mechanism is tha t it is predicted 



to spa wn massive fragments (jWhitworth et al 
1994bh . and hence it affords the possibility of 



forming massive stars sequentially. 

In contrast, RDI dKessel-Devnet fc Burkert 



20031: ISandford et"all Il982l: iBisbas et all \20m 



triggers star formation by compressing pre-existing 
cold, but otherwise gravitationally stable, molecu- 
lar cloud cores. P revio us observational and theo- 
retical studies (e.g. lGritschneder et al. 2009. 2010l) 
suggest that RDI leads to star formation in the 
tips of pillar-like structures, for exam ple as seen in 
the Ea gle nebula (.White et al.J999iV [ Bisbas et all 



lave published a detailed numerical study 
of RDI in initially stable Bonnor-Ebert spheres. 

The nearby HII region RCW 120 is one of 
the most studied HII regions in the Galac- 
tic plane. As observed with Spitzer at 8 fim 
(jChurchwell et all l2006l) . it appears to be an 
almost perfectly roun d bubble with a wel l de- 
fined ionization f r ont. Zavagno et al. ( 2007t ) 
Deharveng et al.l 



( 20091) 



and 



have combined Spitzer 
and 2MASS data with observations at 870/im and 
1.2mm, to analyze the star formation associated 
with RCW 120. They infer that this star for- 
mation has been triggered by the C&C mecha- 
nism. Similar conclusion s have been drawn for 
other regions, e. g. Sh 104 (IDeharveng et al.ll2003l) 



and RCW 79 (jZavagno et al.l 120061 ) where the 
C&C scenario is preferred due to the existence of 
very massive cold clumps. Performing radiation- 
magnetohydrodyna mic simulations of HII regions. 



Arthur et al.l (|201ll ) point out that bubble struc- 
tures like RCW 120 should be quite common. 
However, their simulations do not include self- 
gravity and therefore conclusions on shell frag- 
mentation, clump formation and triggered star 
formation could not be drawn. 

In this paper we show that a distribution of 
massive clumps around an HII region can be ex- 
plained without invoking the C&C mechanism, 
if the HII region expands into a non-uniform 
medium. This scenario does not require a shell to 
be collected and then fragment. Indeed, we find 
that the C&C mechanism is only viable if the cloud 



into which the HII region expands is extremely 
uniform, which seems very unlikely. Observed 
clumpy, shell-like structures like that seen in RCW 
120 are most likely the consequence of pre-existing 
structures in the natal molecular cloud. 

The plan of this paper is the following. In sec- 
tion |2| we describe the algorithm used to generate 
initial fractal molecular clouds, and the numer- 
ical method used to evolve them, including the 
treatment of ionizing radiation. In section |3| we 
describe the resulting HII regions and the mod- 
eled synthetic 870/im observations. We discuss the 
shell and clump masses inferred from the synthetic 
observations in section |4l and compare them with 
the true masses. Our main conclusions are sum- 
marized in section [5] 

2. Initial conditions &: numerical method 

2.1. The generation of a fractal molecular 
cloud 

The initial three-dimensional fractal density 
structure is constructed using an FFT-based al- 
gorithm. The algorithm has three main input pa- 
rameters, (i) the 3D power spectral index n, where 
P{k) oc A;", (ii) the random seed TZ used to gen- 
erate a particular cloud realisation, and (iii) the 
density scaling constant po (see below). We pop- 
ulate the integer modes k = 1,..,128 along each 
Cartesian axis (x,y,z), where k = 1 corresponds 
to the linear box size in one dimension. According 
to Stutzki et al. 1 (|l998i) . the fractal box-coverage 
dimension, D, of a fractal structure embedded in 
three-dimensional space is related to the power 
spectrum by 



P = 3 



(1) 



(see also iFederrath et al.l l2009l ). Thus, defin- 



ing 2? is equivalent to defining the power spec- 
tral index n. Here, we choose setups with P = 
2.4, in agreement with obs ervations of molecular 



clouds in the Milky Way (F algarone et al.l 11991 



Vogelaar & Wakker 



2004; Sanchez et al 



]_994; Stu tzki et al.lll998l : iLee 
2005) . This corresponds to 



n = 3.2. 

After constructing the density fluctuations in 
Fourier space and applying the EFT to give 
p^^^{x,y, z) on a 128'^ grid, the resulting field 
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is scaled using an exponential: 



exp 



(2) 



where po is an arbitrary scaling constant. The re- 
sulting density field has log-normal density proba- 
bility density function (PDF) and a clump mass 
distribution in agreement with observation s (as 
described in IShadmehri &: Elmegree 3 l2Qllh . In 
particular, for a given spectrum of density fluctu- 
ations, changing po allows us to adjust the width 
of the density PDF, i.e. the variance a of the log- 
normal distribution, whilst leaving the underlying 
topology of the density field unchanged. In Figure 
[1] we show the resulting density PDFs for setups 
having equal fractal dimension V — 2A and the 
same random seed TZ, but two different values of 
Po- 

2.2. Initial conditions 

We choose a cloud with total mass of M^^ = 
10"* M(7), and radius of i?„„ 



This re- 



^^g cm 



= 5.0 pc. 

suits in a mean density of p = 5.42 x 10" 
or equivalently n = 1380 cm~^ for molecular gas 
having mean molecular weight /i = 2.35. In this 
paper, we discuss two simulations, both of which 
result in a shell-like structure very similar to the 
one observed in RCW 120. Apart from the scal- 
ing parameter po) the initial clouds are identical, 
i.e. their density fields have the same topologi- 
cal structure. The simulation with po = 1-5 is 
called Run 1, and the simulation with po = 1.0 is 
called Run 2. By fitting their density PDFs using a 
X-squared minimisation method, we estimate the 
corresponding logarithmic density variances, viz. 
CTi = 0.88 {Run 1), and ctz = 1-31 {Run 2). 

We can relate a to the conditions produced 
in turbulent gas by introducing the scaling re- 
lation between the width of the density PDF 
and the tur b ulent Mach number derived by 
Padoan et all (Il99/t) . IPadoan fc Nordlundl (|2002l ) 
and lFederrath et all (|2008h . 



In 1 + b^M 



2^,^2^ 



(3) 



with b ~ 0.5. According to this relation, the den- 
sity field in Run 1 corresponds to a molecular cloud 
with Mach 2.2 turbulence; and Run 2 corresponds 
to Mach 4.3 turbulence. 



Before populating the computational box with 
SPH particles, we shift the point of maximum den- 
sity to the center of the computational domain, 
and position the ionizing star there. Then we cut 
out a sphere with radius i?^^ > centered on the ion- 
izing star. Finally, we partition the computational 
box with a 128^ grid, compute the mass in each 
cell of the grid, and apportion each cell the cor- 
responding number of SPH particles, distributed 
randomly within the cell. 



D = 2,4, p„=1.0 
D = 2,4 
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Fig. 1. — Density PDFs derived for fractal cloud 
having the same fractal dimension D = 2.4 and 
random seed TZ, but different density scaling pa- 
rameters Po = 1.0 and po — 1.5 (black lines). Both 
distributions have been fitted with a log-normal 
(red lines). 



2.3. Numerical method: SPH 



We use the SPH code Seren (jHubber et al 
2011al ). which is well-tested and has already 
been appl ied to many prob l ems in star forma- 



tion (e.g. IWalch et al. 120111: iBisbas et al.l 12011 



Stamatellos et al.ll201l[ ). The ionizing radiation 



is treated with a HEALPix-based adaptive ray- 
splitting algorithm, which allows for optimal res- 
olution of the ionization front in h igh resolution 
simulations (see Bisbas et al. 20091). We emplo y 



the standard SPH algorithm ( Monaghan 19921) . 
The SPH equations of motion are solved with a 
second-order Leapfrog integrator, in conjunction 
with a block time-stepping scheme. Gravitational 
forces are calculated using an octal spatial decom- 
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position tree ( Barnes fc Hutlll986l ). with monopole 
and quadrupole terms and a Gadget- style opening- 



N,. 



^O48.o4±o.i g-1 (jzavagno et al.ll2007[) . 



angle criterion ( Springel et al. 200ll ). We use the 



standard artificial viscosity p rescription 
( Monaghan fc Gin gold' '1983*). moderated with a 



Balsara switch ([Ba lsara 1995). 

The temperature of ionized gas particles is set 
to 10,000 K. The temperature of neutral gas is 
given by a barotropic equation of state, 



(4) 



30 K, 



where the isothermal temperature is 
the critical density is p^rit = 10"^'^ g cm~'^, and 
the adiabatic index is 7 = 5/3. The use of 
^MiN = 30K might influence the fragmentation 
properties of the forming shell. Since the shell 
becomes very dense and should therefore be al- 
lowed to cool further, it might fragment more effi- 
ciently than currently seen in our simulations. For 
this reason will explore a more complicated cool- 
ing function in a future paper. 

We introduce sinks at density peaks above 
PsiNK — 10~^^gcm~'^, provided that the density 
peak in question is at the bo ttom of its local grav- 
itational potential well (see Hubber et al.llioilbl 
in prep.). Since Psm^ 3> Pcrit) ^ condensation 
that is translated into a sink is always well into 
its Kelvin-Helmholtz contraction phase. Once 
formed, a sink is able to accrete gas smoothly 



Hubber et al1l2011bl 



from its surroundings (see 
in prep.) and thereby grow in mass. Sinks are 
identified as protostars, and their radiation can be 
included in the radiative transfer models produced 
with Radmc-SD in the post-processing step (see 
section [331) . 



3. Results 

3.1. SPH simulations 

In this subsection we discuss the results of the 
SPH simulations. In both Run 1 and Run 2 
a source of ionizing radiation emitting A^Lyc ~ 
10^^ s~^ is placed at the center of the cloud. This 
corresponds approximately to an 07.5 ZAMS star 
with a mass of ~ 25M^ and a su rface tempera ture 
of - 40, 000 K (jOsterbrock fc Ferlandl liooj Ta- 
ble 2.3). Thus, the used central source is slightly 
stronger than the central star within RCW 120, 
which is estimated to be a 08.5 - 09 star emitting 



In both simulations we use 2.5 x 10^ SPH par- 
ticles. Thus, each particle has a mass m ~ 
4.0 X 10^'^ Mq and the minimum resolvable mass 
is ~ 0.4 M0. 

In FigurelH we show the initial and final column 
density distributions for the simulations. Note 
that Run 2 has the broader density PDF, and 
hence the more pronounced density contrasts in 
the initial conditions. Both setups develop an HII 
region with a diameter of ~ 5 pc. We compare 
the two simulations at the time when a total mass 
of ^ 500 Mq has been converted into stars, which 
is t^^^ = 0.98 Myr for Run 1 and t^^^ = 0.68 
Myr for Run 2. The black dots in the evolved 
HII regions mark sinks, i.e. protostars. At t^No 
there are 79 sinks in Run 1 and 38 sinks in Run 
2. We stress that the simulations are highly dy- 
namical. New protostars are constantly formed 
and existing protostars continue to accrete at dif- 
ferent rates. In Figure [Sj we plot the mass accre- 
tion rate onto each protostar as a function of its 
mass. The accretion rates of protostars that have 
essentially stopped accreting are arbitrarily set to 
Mmin = 10"^° Moyr"^ For all other sinks there 
is no clear correlation between sink mass and mass 
accretion rate, and we find a mean accretion rate 
of about 10~^ MQyr~^ in both simulations. Run 
2 forms more massive stars than Run 1. 

3.2. Overall bubble structure 

As with RCW 120, we find that the shells 
formed in our simulations are not perfectly spher- 
ically symmetric, but elongated and perforated. 
For example, in the simulations illustrated in Fig. 
[2] the initial cloud has a region of reduced column 
density in the northwest corner, where the ionized 
gas is able to stream out of the HII region and 
open a hole in the shell. 

iDeharveng et al.l ( 2009t ) note the presence of 
radially oriented, filamentary streams below the 
dense, southern part of the shell in RCW 120, and 
suggest that these filaments are due to the pres- 
sure exerted by the ionized gas. Our simulations 
suggest an alternative mechanism for the forma- 
tion of such filaments. As shown in the bottom 



^Both runs could be followed further, but become extremely 
time-consuming and slow once the dense shell is collapsing 
in many places. 
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Fig. 2. — Upper panels: Initial column density distribution in [Mq pc^^] for two representative fractal clouds. 
Both initial conditions were derived with a fractal dimension of 2? = 2.4 and the same random seed TZ, but 
different scaling densities po = 1-5 (Run 1; left column) and po = 1.0 (Run 2; right column). Lower panels: 
The column density distributions at ^^nd after an HII region has formed and part of the molecular cloud has 
been swept up into a shell. The black dots mark sink particles, i.e. forming protostars. 
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Fig. 3. — Mass accretion rate vs. current sink 
mass for Run 1 and Run 2 at t„^,^ . 
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Table 1: Properties of the protostars embedded in 
the shell and the main clumps. Column 1 gives 
the ID of the clump. Columns 2 and 3 give the 
total mass and total number of protostars embed- 
ded within the clumps at <jj,nd • Note that the sink 
masses do not add up to the total mass in proto- 
stars at ijjjND 1 because some protostars are located 
outside the main clumps. 



left panel of Fig. [2l radial streams outside the 
shell, which appear similar to those found in RCW 
120, are developed in Run 1 after ~ 1 Myr. These 
streams are due to the overall gravitational col- 
lapse of the molecular cloud; given a non-uniform 
initial density distributions, the formation of fil- 
aments during gra vitational collapse is inevitable 



1^ K 

(jWalch et al.ll2011l ). We therefore suggest that the 
observed radial filaments could have formed as a 
result of gravitational collapse, independent of ion- 
ization feedback. 

We are unable to identify pronounced pillar- 
like structures at the resolution (19.2") of our syn- 
thetic 870 /im images. The SPH mass distribution 
does show some small pillars and EGGs close to 
the northwest hole, but in general pillars are not 
a prominent feature of HII regions expanding into 
clouds with fractal dimension V = 2.4, because 
- for this fractal dimension - large-scale density 
fluctuations dominate the cloud structure. 

3.3. Synthetic observations 

In a post-processing step we map the SPH den- 
sity distribution onto a three-dimensional grid, us- 
ing kernel-weighted interpolation. We choose the 
grid resolution to be equal to the sp atial reso- 
lution of 19.2" (|Deharveng et al.ll2009l ). achieved 
when observing RCW 120 with APEX-LABOCA 
at 870 /xm. Assuming a constant gas-to-dust ra- 
tio of 100, this immediately gives the dust density 
distribution from the SPH simulation. The only 
modification to this rule is that we assume dust 
to be destroyed at gas temperatures higher than 
1200 kH and therefore the dust density is set to 
a small offset value in these cells; the offset value 
is 10""^ times the minimum dust density in the 
rest of the computational domain. We only take 
into account silicate dust grains, on the assump- 
tion that these dominate the opacity at 870 /im. 
We use t he standard opacity table for this species 
given by iDraine fc Led (|l984l) . 

Using the gridded density distribution, we 
model the transport of continuum radiation 
against dust opacity using RADMC-3E0 (ver- 



^Temperatures higher than 1200 K are exclusively found 
within the HII region. Nevertheless, we tested our results 
against the dust destruction temperature. Increasing it to 
1600 K did not show any effect. 

^http:/ /www. ita.uni-heidelberg.de/ duUemond/software/radmc- 
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sion 0.25; Dullemond, in prep.). In the first 
step RADMC-3D performs a thermal Monte Carlo 
(MC) radiative transfer simulation to determine 
the equilibrium dust temperature distribution. 
The MC method is based on the method of 



Biorkman fc WoodI (|2001l ). but includes various 



improvements, for instance the continuous absorp- 
tion method of Lucy ( 19991 ). The total luminosity 
of all sources within the computational domain 
is distributed amongst Np^^^^ = 10 x 10^ photon 
packages, and N^cat = 5 x 10^ photon packages 
are used to compute scattering events. In one set 
of radiative transfer calculations we only invoke ra- 
diation from the central ionizing star; henceforth 
we refer to these calculations as ionizing source 
only. In addition, we perform radiative transfer 
calculations in which radiation from the newly- 
formed protostars is also includedEI; henceforth 
we refer to these calculations as secondary sources 
included., and distinguish quantities derived from 
these calculations with a superscript 

In the second step, RADMC-3D computes 
isophotal maps at 870 /im, using ray tracing. Fig- 
ure |4] shows the resulting isophotal maps of our 
simulations, as seen in their {x, y)-projection. The 
outer white contour marks the O.lJy/beam cut- 
off value, which we use to define the total mass 
of the shell. The inner white contour marks 
the 0.5 Jy/beam value, which we use to define 
the masses of the clumps. With radiation from 
the ionizing source only (top row of Fig. U), 
the total flux an observer at 1.34 kpc would re- 
ceive at 870 Aim is F*°' = 312 Jy for Run 1 and 
ptot _ J jTQj. ^yj^ 2, If radiation from the sec- 

870 ^ 

ondary sources is included (bottom row of Fig. |4]) 
the resulting total fluxes are signiflcantly higher: 
= 500 Jy for Run 1 and = 760 Jy for 

870 870 

Run 2. The fluxes are higher in the second case, 
because most of the newly formed protostars are 
located in or near the clumps. 

Qualitatively, the synthetic isophotal maps of 
our simulations are very similar to the 870 /im ob- 
servations of RCW 120. In the following, we com- 
pare synthetic and observed images in greater de- 
tail in order to assess the SPH simulation, but also 



3d 

*The luminosity of the sink sources is calculated according 
to their mass, assuming that they are main-sequence stars; 
this probably results in a significant underestimate of their 
luminosities. 



to evaluate the uncertainty of the mass estimates 
for clumps in the shell of RCW 120. 

We use the 0.5 Jy/beam contour to define 
clumps within the shell. In particular, we identify 
the 3 most massive clumps, which we label CI, 
C2, C3, and focus our analysis on their proper- 
ties. The total mass and number of sinks located 
within each clump are listed in Table [TJ The mass 
in gas and dust in the shells and in the individual 
clumps are listed in Table 2, which we describe 
more thoroughly in the next section. 

4. Discussion 

4.1. Clump masses 

In order to compare the simulations with ob- 
servations o f RCW 120, we a d opt th e same tech- 



nique as in Deharveng et al. ( 20091 ) to calculate 



the masses of the shells and of the individual 
clumps. 



IVL. 



100- 



F Z)2 

^ 870 ^ 



5^870 (-^DUST ) 



(5) 



1.8cm2g-i 



Here D is the distance of the source {D 
1.34 kpc for RCW 120), 

(|Ossenkopf &: Hennina Il994l ) is the dust opacity 
per unit mass at 870 /xm, and B^.^^,{T-D\]aT) is the 
Planck function at 870 micron for dust tempera- 
ture Tdust- We adopt Tj^us^ = 30 K, since we do 
not allow the gas to cool below this temperature 
in the simulations. A constant gas-to-dust ratio 
of 100 has been assumed. 

In Figure [5] we show the mass distribution 
derived from the isophotal maps using Eq. [5l 
For reference we overplot the 0.1 Jy/beam and 
0.5 Jy/beam contours from Fig. 21 which define 
- respectively - the shell and the most massive 
clumps. For Run 1 we calculate a total shell mass 
of 1818 Mq with radiation from the ionizing source 
only, and 2928 Mq with radiation from secondary 
sources included. Both of these estimates are sim- 
ilar to the shell mass of 1, 100 M© found in RCW 
120 with T^^jsT = 30 K. For Run 2 the correspond- 
ing shell masses are 3013 M© and 4420 M©. 

4.2. Massive clumps without C&C 

We divide the shells up into clumps using the 
0.5 Jy/beam contour on the synthetic isophotal 
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870 micron emission 870 micron emission 




X [pc] 




Fig. 4. — 870 micron emission calculated using RADMC-3D. The left column shows the images for Run 1, 
and the right column for Run 2. The top row shows the images obtained with radiation from the ionizing 
source only, whereas the bottom row shows the images obtained with radiation from secondary sources 
included. The white contours are set to fluxes of 0.1 Jy/beam and 0.5 Jy/beam in all images. 



8 



Mass distribution 



Mass distribution 



0.6 



0.2 



^ 






2 
[pc] 



1.5 



^ 




Mass distribution 



Ivlass distribution 



0.5 



0.0 



^ 



4 





2 
[pc] 



0.5 



0.0 



^ 



2 
[pc] 



Fig. 5. — Mass distribution derived from tlie tliermal dust emission at 870 /im shown in Fig. |4l For 
comparison we overlay contours marking 0.1 Jy/beam and 0.5 Jy/beam from Fig. |4l As in Fig. |4j tlie top 
row lias been calculated with radiation from the ionizing source only, and the bottom row with radiation 
from secondary sources included. The three largest clumps, which are discussed in the text, are marked with 
CI, C2, and C3. 



maps. These clumps typically have masses be- 
tween a few and a few hundred M-q, depending 
on the initial cloud density and density PDF. In 
Fig. [S]we identify the three main clumps formed 
within each shell and list their properties in Ta- 
ble [2] All these clumps are sufficiently massive to 
possibly spawn massive protostars. 

The clumps are nicely aligned with the form- 
ing dense shell and the ionization front. However, 
our simulations clearly show that the formation 
of separate clumps in these simulations is not due 
to C&C. At no time do we observe the formation 
of a coherent shell, which grows to become grav- 
itationally unstable and then undergoes fragmen- 
tation. The formation of clumps simply reflects 
the dumpiness of the underlying initial molecu- 
lar cloud structure. Density variations are nat- 
urally enhanced by the ionizing radiation, which 
penetrates low-density regions much more easily 
than high-density regions. The consequent heat- 
ing of the low-density regions increases their pres- 
sure, and this acts both to disperse the low-density 
gas, and to compress the remaining high-density 
clumps. The main clumps in Fig. [5] can clearly be 
related to structures in the initial column density 
profile in Fig. [5] 

The masses estimated for the three main 
clumps are strongly dependent on whether the 
radiative transfer modeling includes the radiation 
from secondary sources or not (see Tableland sec- 
tion l4.3p . The variations can be as high as a factor 
of 8, as is the case for C3 in Run 1, for which the 
estimated mass increases from 61 M0 to 493 M© 
when the radiation from secondary sources is in- 
cluded. This is because the extra heating from 
newly-formed protostars makes the dust in their 
vicinity hotter, and therefore more material falls 
within the 0.5 Jy /beam threshold. 

4.3. How reliable is the mass distribution 
obtained from 870 /im fluxes? 

In this section we compare how well the actual 
mass distribution in the simulations is recovered 
by applying Eq. [5] to the synthetic isophotal maps 
at 870 /im. This comparison can provide useful in- 
sights into the reliability of clump mass estimates 
from observational data. 

We define four masses, which are listed in Ta- 
ble El Mg,o is the mass obtained using Eq. [5] 



on synthetic isophotal maps calculated with ra- 
diation from the ionizing source only, and M* 
is the mass obtained using Eq. [S] on synthetic 
870 /im isophotal maps calculated with secondary 
sources included. Likewise, M^j^^b is the actual 
mass falling within a shell or clump on synthetic 
870 /im isophotal maps calculated with radiation 
from the ionizing source only, whilst M* is the 

o J 7 TRUE 

corresponding quantity when radiation from sec- 
ondary sources is included in the radiation trans- 
fer modeling. M^j^^^, and M*^^^ are obtained by 
integrating the surface density of SPH particles 
over the area covered by the shell or clump (i.e. 
the area inside the 0.1 Jy/beam and 0.5 Jy/beam 
contours, respectively). 

The resulting fractional errors. 



1-^870 -^TRUE I 
M 



and 



I 870 



M* 

TRUE I 



(6) 



(7) 



are shown in Fig. [6l The shell masses obtained 
using Eq. [5] on synthetic 870 /im isophotal maps 
are always lower than the actual mass, typically 
by a factor ~ 2, irrespective of whether the radi- 
ation from secondary sources is included or not. 
The masses of individual clumps obtained using 
Eq. [S] on synthetic 870 iim isophotal maps are 
sometimes higher and sometimes lower than the 
true mass; for massive clumps (> 100 M©) the 
fractional error is always less than two, but for 
lower-mass clumps it can be very large, and such 
masses should probably not be trusted. 

In general, masses obtained using Eq. [5] on 
synthetic 870 isophotal maps calculated with 
radiation from secondary sources are more accu- 
rate. However, we note that the inclusion of ra- 
diation from secondary sources also increases the 
dust temperature, and this should be taken into 
account when applying Eq. [5l In particular, for 
the three clumps identified in Run 2, Radmc-3d 
gives dust temperatures of about 40 K. This causes 
us to overestimate the clump masses M*^^ by a 
factor of 



^.™(30K) 



1.4. 



(8) 



10 



10.0 



ID 


Mq 


M* 

870 

Mo 




TRUE 

Mq 






Run 


1 




Shell 


1818 


2928 


3031 


4450 


CI 


33 


91 


5.0 


23. 


C2 


75 


146 


114 


262 
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61 
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336 






Run 


2 




Shell 


3013 


4420 


6291 


6853 


CI 


193 


302 


129 


228 


C2 


74 


184 


114 


270 


C3 


160 
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89 
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Table 2: Column 1 identifies the different struc- 
tures analysed. Columns 2 and 3 give mass esti- 
mates derived from the synthetic 870 /xm images 
which were calculated with radiation from the ion- 
izing source only, M^^„, and with radiation from 
secondary sources included, M*^^ (see Fig. 5). 
Columns 4 and 5 give the true masses, obtained 
by integrating the SPH column density over the 
area inside the O.lJy/beam contour for shells, and 
the 0.5Jy/beam contour for the clumps. M^-^-^J-^ 
is the mass within these contour levels calculated 
with radiation from the ionizing star only, whereas 
M* is the mass within these contour levels cal- 

TRUE 

culated with radiation from the secondary sources 
included. Both, A/true ^"^^ -^^toue ^^ccount only 
for the mass of gas and dust and do not include 
embedded protostars. The mass in sinks is listed 
in Table 1. 
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Fig. 6. — Fractional error in the masses estimated 
for shells and clumps obtained using Eq. [S]on syn- 
thetic isophotal maps with and without radiation 
from newly- formed protostars (see Eqs. |6]and[7l). 



5. Conclusions 

We have performed high resolution SPH sim- 
ulations of HII regions expanding into fractal 
molecular clouds, and compared synthetic 870 /im 
isophotal maps of these simulations (obtained us- 
ing RADMC-3D) to 870 /xm observations of the 
well-studied, galactic HII region RCW 120. 

We conclude that finding massive clumps and 
sites of high mass star formation within an ex- 
panding shell formed around an HII region pro- 
vides neither proof nor need for the C&C sce- 
nario. We suggest that a more accurate descrip- 
tion of the triggering mechanism is given by the 
Enhancement of initial Density substructure and 
simultaneous Global Implosion (EDGI) of clumps 
due to the more efhcient heating and dispersal of 
low-density regions in between denser structures. 

The intrinsic dumpiness of the forming dense 
shell reflects the pre-existing cloud density struc- 
ture. To first order, regions which later form 
clumps can already be identified in the initial con- 
ditions as areas of enhanced column density. 

Overall, we find good agreement (to within a 
factor of 2) between the actual mass distribution 
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and the mass distribution inferred from thermal 
dust emission, for shells and clumps having masses 
greater than 100 M©. In particular, masses ob- 
tained using Eq. [5] on synthetic 870 ^m isopho- 
tal maps always underpredict the overall mass of 
the the large-scale shell structure, whereas for in- 
dividual clumps both over- and under-estimates 
are possible, and usually the result is accurate to 
within a factor of two. 

We note that filamentary streams that arc ori- 
ented towards the center of the Hll region and are 
located outside of the ionization front arc prob- 
ably not related to, or caused by the HII region 
itself. We find that such streams form in our sim- 
ulations due to the overall gravitational collapse of 
the cloud. Consequently, these dense streams tend 
to form behind regions where most of the ionizing 
radiation is trapped by a dense shell of gas and 
dust. The presence of filamentary streams may 
hint at the age of the cloud. 
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